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ABSTRACT 

We present and discuss measurements of the gas-phase metalhcity gradient in gravitationally lensed 
galaxies at z = 2.0 — 2.4 based on adaptive optics-assisted imaging spectroscopy with the Keck II 
telescope. Through deep exposures we have secured high signal to noise data for four galaxies with 
well-understood kinematic properties. Three galaxies with well-ordered rotation reveal metalhcity 
gradients in the sense of having lower gas-phase metallicities at larger galactocentric radii. Two of 
these display gradients much steeper than found locally, while a third has one similar to that seen in 
local disk galaxies. The fourth galaxy exhibits complex kinematics indicative of an ongoing merger 
and reveals an "inverted" gradient with lower metalhcity in the central regions. By comparing our 
sample to similar data in the literature for lower redshift galaxies, we determine that, on average, 
metalhcity gradients must flatten by a factor of 2.6 ±0.9 between z = 2.2 and the present epoch. This 
factor is in rough agreement with the size growth of massive galaxies suggesting that inside-out growth 
can account for the evolution of metalhcity gradients. Since the addition of our new data provides 
the first indication of a coherent picture of this evolution, we develop a simple model of chemical 
evolution to explain the collective data. We find that metalhcity gradients and their evolution can be 
explained by the inward radial migration of gas together with a radial variation in the mass loading 
factor governing the ratio of outfiowing gas to the local star formation rate. Average mass loading 
factors of < 2 are inferred from our model in good agreement with direct measurements of outflowing 
gas in z ~ 2 galaxies. 

Subject headings: galaxies: high redshift - galaxies: evolution 



1. INTRODUCTION 

Metalhcity is one of the most fundamental properties 
of a galaxy. Gas-phase metalhcity (hereafter referred 
to simply as metalhcity) is governed by the cumulative 
history of baryonic assembly: gas accretion, star forma- 
tion, and gas outflow. Metalhcity is therefore a tracer of 
the processes responsible for galaxy formation and evo- 
lution. The vast database of information made available 
by the Sloan Digitized Sky Survey (SDSS) has revealed 
that metalhcity is correlated with galaxy stellar mass and 
anticorrel ated with star formati on rate ([T remont i et al.l 
[2004: M annucci et al . 2010; Lar a-Lopez et a l. 201Q1. The 
dependence on stellar mass is explained as the effect of 
gravity: more massive galaxies lose a smaller fraction of 
their metal-enriched gas in outflows and are more eas- 
ily able to re-accrete lost metals due to their stronger 
gravitational potential (Tremonti ct al. 2004). The de- 
pendence on star formation rate is primarily an effect of 
accretion: increased accretion of metal-poor gas drives 
higher star formation rates and lowers the overall metal- 
hcity (e.g. iDave et al.ll20lH ). 

The spatial distribution of metalhcity and its evolu- 
tion with time is sensitive to the assembly history of 
galaxies. Careful measurements in the local universe 
show that all disk galaxies exhibit negative radial metal- 
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licity gradients, with lower metalhcity at larger galac- 
tocentric radi us (e .g. Vila-Costas & Edmunds_ 1992t 

'va n Zee et allll998t ICMTsidere et al.H2000l: (Rupteet alj 
I2010bl) . In many cases the gradient flattens at large 
radius indicating efficient radial mixing, possibly in- 
duced by interactions or cycling between t he disk and 
halo i n a "galactic founta in" process (e.g. iWerk et aD 
120111: iBresolin et~aLl I2012D . More attention has been 
given recently to measurements of the time evolution 
of metal gradients and the implications for galaxy for- 
mation. Models of inside-out disk growth predict that 
gradients shoul d be initially steep and become flatter at 
later t imes (e.g . Prantzos & Boissicr 2000^'Magrini et all 
120071: iFu et alj 2009; Marcon-Uchida ct al. 2010). In 
contrast, numerical simulations show that merger-driven 
growth will rapidly flatten existing meta hicity gradi- 
ents, w hich then become steeper with time (jRupke et al.l 
l2010al ). This is supported by nearby observations which 
reveal flatter gradients in inte racting galaxies com pared 
to an isolated control sample (jRupke et al.ll2010b[) . 

Several groups have now begun to consider the time 
evolution of radial metalhcity gradients through obser- 
vations. One approach is to examine the situat ion in 
detail within the Milky Way. IMaciel et all (|2003D have 
inferred that the metalhcity gradient must have been 
steeper in the past by analyzing the properties of stel- 
lar systems with different characteristic ages. However, 
this method is only feasible in a small number of local 
galaxies. The alternative approach, adopted here, is to 
measure the gradient in resolved data for galaxies seen at 
high redshift and to then compare these measures with 
those at lower redshift. Although a challenging task, this 
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was first attempted in a previous paper where we showed 
the metallicity gradient in a gravitationally lensed galaxy 
at z = 2 was sign ificantly steeper than in local disk 
galaxies (jJones et a l. 2010b). Another lensed g alaxy at 
z = 1.5 yielded similar results (jYuan et al.ll201lT ). further 
suggesting an inside-out mode of galaxy growth. How- 
ever, some studies at high redshift have suggested that 
many massive galaxies may have "inverted" (positive) 
gradients wit h low er metallicity in the central regions 
(|Cresci et alllMoHOuevrel et aLllMl . This would im- 
ply a radically different mode of growth. One issue with 
this approach is the technical difficulty of measuring a re- 
liable gradient in seeing-limited data which offers a rel- 
atively poor spatial resolution of ~ 4 kpc FWHM, ap- 
proximately twice the typical half-light rad ius of an 
galaxy at z = 2 — 3 (jBouwens et al.l [20041 ) . Such mea- 
surements cannot provide good spatial sampling except 
for the very largest sources. Clearly it is advantageous to 
use adaptive optics as well as the angular magnification 
afforded by gravitational lensing. 

A further issue concerns the fact that some high red- 
shift metal gradients reported in the literature are mea- 
sured from a limited set of emission lines which may 
be affected by shocks and/or AGN (e.g. [Nii]/Ha, 
[Oiii]/H/3). Shocks can easily produce a signature indica- 
tive of an inver ted metallicity gradient in the [Nil] /Ha 
hne ratio (e.g. iWestmoauette et"ani2012[ ). while AGN 
can mimic an inverted gradient in [O ili]/H/3 and a stron g 
negative gradient in [Nil] /Ha (e.g. iWright et al.ll2010l ). 
Hence multiple emission line diagnostics are required to 
distinguish variations in metallicity from the effects of 
shocks and AGN fo llowing methods such as those of 
[Baldwin et al.l ()1981t ). While observations of multiple 
emission lines required for such diagnostics are challeng- 
ing to obtain, they are necessary for robust metallicity 
gradient measurements. 

This paper is concerned with increasing the number 
of galaxies at high redshift with reliable high-resolution 
metallicity measurements. Gravitationally lensed galax- 
ies are ideally suited to this purpose, as the strong magni- 
fication provides bot h an inc reased flux and higher source 
plane resolution (e.g. I Jones et al. 2010a.b.'). In this paper 
we present spatially resolved metallicity measurements of 
four lensed gala xies at z = 2 — 2.4 i ncluding one previ- 
ously studied bv I Jones et al.l ()2010b[ ). For the first time 
this permits us to consider the overall evolution of the 
metallicity gradient by utilizing lower redshift data. In §2 
we describe the source selection, gravitational lens mod- 
els, observations, and data reduction. In §3 we discuss 
the kinematic properties of each source which provides 
valuable insight into variations we see within our sample. 
In §4 we present resolved metallicity measurements. In 
§5 we construct and discuss a simple chemical evolution 
model to describe the origin of metallicity gradients in 
the context of the growing body of data we assemble for 
comparison purposes. The details of this model are given 
in an Appendix. In §6 we discuss the results in the con- 
text of our chemical evolution model. Finally, in §7 we 
summarize the main conclusions and implications of the 
results. 

Throughout this paper we adopt a ACDM cosmology 
with Ho= 70 km s^^ Mpc~\ r2Af=0.30 and f]A=0.70. 
At z=2.2, 0.1 arcsec corresponds to 830 pc and the age 
of the universe was 2.9 Gyr. All magnitudes are in the 



AB system. 

2. OBSERVATIONS AND DATA ANALYSIS 
2.1. Source Selection 

Our sources were selected as ideal targets for resolved 
metallicity measurement based on their redshift, nebu- 
lar emission line intensity, magnification, and the pres- 
ence of a suitable star for tip/tilt correction. Our lensed 
sources were selected to lie at z = 2.0 — 2.4 such that 
the diagnostic emission lines Ha, H/3, [Nil], and [Oiii] 
can be observed in the H and K atmospheric windows 
where the Keck II adaptive optics (AO) system provides 
the best Strehl ratio. The exception is SDSS J1206 for 
which H/3 lies in a telluric absorption band and is inac- 
cessible with OSIRIS. Prior to AO-assisted observation, 
we secured a near infrared spectrum of each source to en- 
sure that the nebular emission lines are sufficiently bright 
and relatively free o f contamination from OH sky lines 
(|Richard et al.l[20TTh . We require that each gravitational 
lens system has an accurate mass model constrained 
by the positions of multiply-imaged background sources 
with known redshift, essential for the source plane recon- 
struction discussed in §2.21 Finally, we require a suitably 
bright and nearby star to provide tip/tilt correction for 
the adaptive optics system. The observational details of 
each galaxy are given in Table [TJ Below we provide a 
brief overview of the sources described in this paper. 

MACS J0744+3927 (hereafter J0744) is a galaxy clus- 
ter lensin g system di scovered by the MAssive Cluster 
burvey fEbeling et al.ir200lD . The z = 2.21 arc was iden- 
tified in foUowup optical spectroscopy of the cluster and 
we secured the near-infrared spectrum a,s part of our 
screening program of bright ^ > 2 arcs (|Richard et al.l 
1201 Ih . We previously observed the Ha emission line 
with OSIRIS and described the source kinematics in 
iJones et al.l ()2010aD . 

SDSS J1038+4849 (hereafter J1038) is a group-scale 
lensing system discovered in the CASSOWARY survey 
for wide-separation gravi tational lenses in SDSS imag- 
ing (jBelokurov et al.ll2009 ;]) . We obtained a near- infrared 
spectrum of the z = 2.20 arc using NIRSPEC on Keck 
which showed bright nebular emission lines suitable for 
resolved spectroscopy with OSIRIS. We show in Sj3]that 
the source is a merger of at least two systems with stellar 
mass ratio of (6±3):1. Therefore in much of the analysis 
we consider the two components separately. The com- 
ponents are divided into the region with y-axis > — 1 
kpc in Figure ( "North" region) and y-axis < — 1 kpc 
("South"). The northern region is brighter in optical 
WFPC2 imaging and the southern region is brighter in 
IRAC. 

SDSS J1148+1930 (hereafter J1148) is comprised of a 
background source at z = 2.38 lensed into a nearly com- 
plete Einstein ring by a massive galaxy at z = 0.44. The 
system was discovered in the SDSS by iBelokurov et al.l 
(|2007l ) , and n ear-infrared lon g slit sp ectroscopy of the arc 
is reported in Hainline et al.l ([200 9^). 

SDSS JI2O6+5142 (hereafter J1206) is a g roup-scale 
lensing system discovered in the SDSS by iLin et all 
(|2009f) . We previously obtained resolved metallicity mea- 
surements of this source from observations of Ha, [Nil], 
and [Oiii] with OSIRIS (IJones et al.ll20l0bl) . J1206 has 
an unusually steep radial metallicity gradient compared 
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Name 


z 


Coordinates 


Dates 


Filter 


Lines 




FWHM 


FWHM 
















(PSF) 


(source plane) 




J0744 


2.21 


07:44:47.9 +39:27:26 


27 Nov 2008 


Kn2 


Ha, [Nil] 


14400 


O'.'ll 


0.3 X 0.8 kpc 


16 + 3 








20 Feb 2011 


Hn2 


[Oiii], H^ 


10800 


0'.'08 


0.3 X 0.7 kpc 




J1038 


2.20 


10:38:41.8 +48:49:19 


20 Feb 2011, 12 Mar 2011 


Kn2 


Ha, [Nil] 


4800 


0'.'14 


0.4 X 1.6 kpc 


8.4 + 0.7 








20 Feb 2011 


Hn2 


[Oiii], H^ 


1200 


0'.'14 


0.3 X 1.7 kpc 




J1148 


2.38 


11:48:33.3 +19:29:59 


20 Feb 2011 


Kc4 


Ha, [Nil] 


7200 


O'.'ll 


0.6 X 0.9 kpc 


10.3 + 5.0 








20 Feb 2011, 12 Mar 2011 


Hbb 


[Oiii], H^ 


3600 


0'.'08 


0.6 X 0.9 kpc 




J1206 


2.00 


12:06:01.7 +51:42:30 


19 May 2010 


Knl 


Ha, [Nil] 


3600 


0'.'18 


0.5 X 3.0 kpc 


13.1 + 0.7 








19 May 2010 


Hnl 


[Oiii] 


1800 


0'.'33 


0.6 X 3.5 kpc 





TABLE 1: Log of observations. 



to local galaxies, motivating high-resolution studies of 
additional sources presented in this paper. We include 
this source in the analysis and discussion of the present 
paper for completeness. 

2.2. Gravitational Lens Models 

Accurate models of the foreground mass distribution 
are essential for correcting the gravitational lensing dis- 
tortion and reproducing the source-plane morphology 
of background galaxies. Here we briefly describe the 
adopted procedure. We use the Lenstool program 
(jKneib et al.lll993l : lJullo et al"1l2007f ) to parameterize the 
mass distribution of the lens, using the positions of 
multiply-imaged background sources as constraints. The 
mass models used for J0744 and J1 206 are described in 
Limousin et al.(in preparation) and lJones et al.l ()2010bO 
respectively. We constructed lensing mass models for the 
two other systems based on available HST imaging and 
multiple images with spectroscopic redshifts. In the case 
of J1038, w e use the two trip le systems at z = 0.9657 and 
z = 2.198 (|Belokurov et al.l 2009) to constrain a double 
pseudo-isothermal elliptical (dPIE) dark matter profile 
(inclu ding core rcore and cut rcut radii, see lJullo fc KneibI 
l2009l for more details). Cut radii are not constrained by 
strong lensing and so we use a fixed value rcut — 800 
kpc typical of galaxy groups. For J1148 we use the sys- 
tem of 4 images at z = 2.38 forming the 'Horseshoe' as 
constraints f or a dPIE mass d istribution, following the 
modelling of iDve et al.l ()2008f ). In both lens systems, 
individual group galaxies detected in the optical bands 
were added as low-scale dPIE mass components perturb- 
ing the model, a pro c edure s imilar to the one used for 
J1206 in iJones et all ()2010b[ ). The best fit parameters 
of the dPIE potentials for these previously unpublished 
cases are summarized in Table [2] 

Lenstool also determines the transformation from 
image plane to source plane position, and we use this 
mapping to reconstruct all observations onto a uniform 
grid in the source plane. Intrinsic source-plane morpholo- 
gies of each galaxy are shown in Figure [2l Lensing mag- 
nification is calculated as the ratio of Ha fiux in the 
image plane to that in the source plane, or equivalently 
the ratio of area subtended by Ha emission. In Table [T] 
we give the typical magnification and la lens model un- 
certainty for the regions observed with OSIRIS. In the 
case of J1206, multiple images are covered with OSIRIS 
such that the total magnification of images within the 
field of view is ~ 20. 

2.3. Integral Field Spectroscopy 

Spectroscopic observations of the targets in Table [1] 
were taken with the OH-Suppressing Infra-Red Imaging 



Parameter J1038 J 1148 

Xc" -2.4+0.7 -0.11+0.11 

Yc ° -2.6+0.5 0.47+0.13 

e'' 0.76+0.03 0.00+0.02 

r 31+1 134+12 

a'' 576+12 654+78 

rcore" 34.1+2.9 30.9+10.2 

rent " [800] [800] 

" Position from the brightest galaxy of the group (arcseconds). 
EUipticity. 

Position angle (degrees). 
Velocity dispersion (kms~^). 
Core and cut radii (kpc). 

TABLE 2: Best fit parameters for the previously unpub- 
lished lensing potentials (see text). 

Spectrograph (OSIRIS; ILarkin et al.l 120061) m conjunc- 
tion with the laser guide star assisted adaptive optics 
(AO) system (IWizinowich et al.l [20061 ) on the Keck II 
telescope. These were carried out during four separate 
observing runs with seeing ranging from 0'.'7 to I'.'S and 
clear conditions as summarized in Table |TJ A suitably 
bright star {R< 17) within 55 arcseconds of each target 
was used for tip-tilt correction. We used the 100 mil- 
liarcsecond pixel scale in all observations which provided 
a field of view between l'.'6x6'.'4 and 4'.'5x6'.'4 depend- 
ing on the filter. Observations of each target were done 
with an AB observing sequence, dithering by ~ 2 — 3 
arcseconds to keep the target within the integral field 
unit. The fields of view for Ha and [Nil] are shown for 
reference in Figure |TJ Fields of view for [Oiii] and H/3 
are essentially identical except for J1038, where only the 
northern region was observed. 

Our data reductio n methods a r e essen tially identical 
to those described in IJones et al.l (|2010b[ ). We used the 
OSIRIS Data Reduction Pipeline (ODRP; ILarkin et all 
l2006f ) to perform dark and bias subtraction, cosmic ray 
rejection, wavelength calibration, and to assemble the 
3D data cubes. S ky subtraction was done with the IDL 
code described in iDaviesI (|2007l ) using temporally adja- 
cent exposures as sky reference frames. Final data cubes 
were combined using a cr-clipped mean. Flux calibration 
for J 1206 and the K-b a nd obs ervati ons of J0744 are de- 
scribed in IJones et al.l (poTOb) and Jone s et al.l (|2010a[ ) 
respectively. For the r emaining ob s ervatio ns, we follow 
the same method as in IJones et al.l (|2010al ) using obser- 
vations of the tip/tilt reference stars to calibrate the ab- 
solute fiux. We additionally check the flux calibration 
with observations of the UKIRT infrared standard star 
FS 26 taken at a different time on the same nights. Flux 
calibrations derived from FS 26 agree to within 15% sug- 
gesting the systematic uncertainty in flux calibration is 
< 15%. 
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SDSSJ1148+1930 



SDSS J1206+5142 
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FIG. 1: Color HST images of the sample. The white box in each image shows the OSIRIS field of view used to observe 
Hq and [Nil] emission fines. R/G/B channels of each image are respectively ACS F814W/F555W/WFC3 F390W for J0744, 
WFPC2 F814W/F606W/F450W for J1038 and J1206, and WFC3 F814W/f'606W/F475W for J1148. 



2.3.1. Extinction and Star Formation Rate 



In this subsection we determine the extinction and star 
formation rate of each galaxy from their integrated spec- 
tra. We determine the total flux of nebular emission lines 
by summing all pixels within 250 kms^^ of the systemic 
redshift and within ~ O'.'S of regions where the arc is 
detected in HST imaging. The results are shown in Ta- 
ble [3l Balmer line ratios Hq;/H/3 are used to determine 
the dust extinction assuming a iCalzetti et al.. (,2000.) red- 



dening curve, such that 



E(B-V) = 1.695 log 



2.86 



(1) 



where Ha/H/? — 2.86 is the intrinsic line ratio for case B 
recombination. The extinction in V band and in Ha are 
then given by 

Av = RvE(B- V) =4.05E(B- V) (2) 

Ana = 3.33 E(B-V) (3) 
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FIG. 2: Source plane reconstruction of the lensed galaxies. Each row shows the morphology of one galaxy in (from left to 
right) rest-UV continuum, Ha emission, and rest-optical continuum. No rest-optical continuum image is available for SDSS 
J1038 so we show the reconstructed IRAC 3.6/xm image. In all sources the ionized gas morphology traced by Ha is similar to 
that of the rest-UV continuum, shown as contours of HST/ACS F814W intensity on the Ha maps. 
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Name 


Ha 


[Nil] 


[Oin]5007 


[OlIl]4959 


H/3 


J0744 


20 ± 2 


9±2 


12 ±2 


6±2 


6 ± 2 


J1038 (North) 


66 ± 5 


5 ± 5 


100 ±5 


42 ±8 


19 ± 7 


J1038 (South) 


71 ±6 


0± 7 








J1148 


67 ±3 


7±3 


41 ±2 


8±2 


10 ±3 


J1206 


236 ±7 


51 ± 7 


230 ±11 


108 ± 14 





TABLE 3: Total observed emission line fluxes. All values 
are in units of 10~^^ ergs~^ cm~^. Uncertainties are la deter- 
mined from noise in the integrated spectra and do not include 
systematic errors in absolute flux calibration, which are typi- 
cally < 15%. For J1038 we report the northern and southern 
regions separately; the southern region was not observed in 
[Oiii] and H^. 

for the ICalzetti et al.l (j2000D reddening curve. Star for- 
mation rates are computed from total Ha flux cor- 
rected for extinction and lensing ma gnificati o n (Ta ble [J), 
with the conversion factor from Ke nnicut^ ()1998f ) for a 
Salpeter IMF. The resulting E(B-V), Ay, SFR, and la 
uncertainties for each source are listed in Table ID We 
note that the lower bound on SFR is required to be 
greater than or equal to the Ha-derived SFR assuming 
no extinction. is not observed for the southern re- 
gion of J1038 and so we give only the extinction-free 
SFR as a lower li mit. For J 1206 we use the Balmer line 
ratios reported by iHainline et al.l ()2009[ ) to derive E(B- 
V) = 0.30 ± 0.12 and apply this extinction correction to 
the Ha flux measured with OSIRIS. 

2.4. Photometry and Stellar Mass 

For an appropriate comparison with lower redshift 
data, it will be helpful to estimate the masses of our 
target galaxies. We have determined stellar masses for 
each galaxy from multi- wavelength photometry and stel- 
lar population synthesis modelling. All of our sources 
have HST imaging in multiple optical filters (e.g. Fig- 
ure H]) as well as Spitzer IRAC channels 1 and 2. All 
galaxies except J1038 additionally have high resolution 
near- infrared imaging from HST WFC3/IR. We deter- 
mine the photometry in matched apertures and then 
scale to the total flux. All data are first smoothed to 
match the point spread function of IRAC channel 1, 
which has the lowest spatial resolution. We then mea- 
sure the flux in each filter within an aperture selected to 
encompass the majority of the arc while avoiding con- 
tamination from nearby sources. The total flux is deter- 
mined from high-resolution optical HST images and all 
aperture fluxes are scaled by the same factor to obtain 
the total flux. Finally, near-infrared fluxes are corrected 
for nebular emission by subtracting the total line fluxes 
(Table [3]). We fit the photometric dat a with stellar pop- 
ulation models using the code FAST (|Kriek et al.ll2009D 
with an exponentially declining star formation history, 
Salpeter initial mass function, and .Calzetti et al. (_2000) 
reddening law. Redshifts are fixed to the spectroscopic 
values and metallicity is allowed to vary within the spec- 
troscopic values for each arc (e.g. Table [S]). We allow a 
range of star formation timescales 10'' — 10^° yr. Ages are 
restricted to be > 50 Myr (approximately equal to the 
dynamical timescale for these sources; see 331) and less 
than the age of the universe (~ 3 Gyr at 2: = 2 — 2.4). 
Dust extinction is allowed to vary within the ilcr range 
of Ay in Table |4] and can be as low as half the allowed 
minimum value, since stellar continuum is often inferred 



to have a factor of c ± 2 lower Ay than Hii regions (e.g. 
iNewman et~ani2012f) . In all cases we require Ay > 0, 
and for the southern region of J1038 we assume Ay < 4. 
Finally, we consider only star formation histories which 
produce a SFR consistent to within la of the Ha-derived 
values. The best fit stellar mass and 68% confidence lev- 
els within these constraints are given in Table SI 

Correcting photometric measurements for contamina- 
tion by nebular emission lines may be important for accu- 
rately determining stellar masses of star forming galaxies. 
Our sources have emission lines with rest-frame equiva- 
lent widths Wo ~ 50 A for [Oiiijsooy and Ha. At red- 
shift z = 2 this corresponds to roughly 10% of the flux 
measured in a broad-band filter. Indeed, [Olll] and H/3 
contribute ~ 20% of total H-band flux measured for the 
lensed sample. To examine the effect on derived stel- 
lar mass, we apply the same stellar population synthesis 
method (and the same constraints) without correcting 
for nebular emission lines. This results in best-fit stellar 
masses which are higher by -1-0.21 dex for J0744, -1-0.002 
dex for J1148, and +0.05 dex for J1206. In general the 
correction is not large and is consistent in all cases within 
the uncertainties, however in the case of J0744 the stel- 
lar mass is overestimated by 60% if nebular emission is 
not accounted for. In summary, ignoring the nebular 
emission contribution would bias the stellar masses to 
somewhat higher values, but not significantly affect any 
results in this paper. 

2.5. Gas Fraction 

The gas mass in star forming galaxies at redshift z ^ 2 
comprises a signific ant fra.ction 50% of thei r tota. 1 bary- 
onic mass (Tacconi eFal] 120101: iDaddi et all [20101) . Ide- 
ally we would like to determine the gas mass from direct 
observations of atomic Hi and molecular CO emission. 
However, in the absence of such data we can estimate 
total gas mass from the Kennicu tt-Schm idt relation. We 
adopt the best-fit relation of KennicuttI (|T998i ) , 

Egas = (4.0 X 10^ Ssfr)°-^' M0PC-2 (4) 

where Ssfr is in units of M0yr^^kpc^^. The baryonic 
gas fraction is then defined as 



We calculate the gas mass and gas fraction of all three 
non-merging galaxies using Ha-derived star formation 
rates and diameters (Table HJ to determine surface den- 
sities. The results are given in Table 21 The inferred 
gas fractions of 0.4 — 0.8 are in goo d agreement with di- 
rect measurement s at high redshift (iTacconi et al.ll20Tol : 
IDaddi et al.ll2010t ). Dynamical masses derived in ^also 
suggest fgas ~ Mdyp-M, _ _ Q g j.^^. ^j^^ a,rcs, in 

excellent agreement. 

In cases where high resolution i/^T infrared imaging is 
available, we calculate the spatially resolved gas fraction 
of the lensed galaxies. For each source-plane pixel we 
determine Ssfa from Ha intensity and the extinction 
factors in Table [H and calculate Egas from Equation [H 
The stellar mass density is calculated for the same pix- 
els from the HST F160W image, which we reconstruct 
in the source plane and smooth to match the resolution 
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J0744 



J 1038 North 



logAf* (Mq) 


^•*-0.i 




E(B-V) 


0.13 ±0.27 


0.17 ±0.32 


Av 


0.53 ± 1.20 


0.67 ± 1.30 


SFR (Mgyr-i) 






Mgas (109 Mq) 






/gas 


42+0-40 




log A/halo (Mq) 


11.6 


11.4 


D (kpc) 


2.4 


13.3 


AV (kms-1) 


252 ± 33 


160 ± 10 


a (kms~^) 


89 ±26 


82 ± 22 


Mayn (lOio Mq) 


1.1 




Ql sini 


0.5 





J1038 South 
TTTT^ 



9.9"* 



> 24 



11.8 



J1148 



9.9 
0.73 ±0.26 
2.94 ± 1.05 

21Q+167 



J1206 



-0.2 



27-^ 



7+14 
-18 

'-0.15 

11.9 
2.7 

148 ±3 
90 ±33 

1.3 

0.3 



10.1] 
0.30 ±0.12 
1.22 ±0.47 

68^ 
14^ 



J +44 
-24 
1+3 
-4 



"S2+''-25 

'J-s^-o.ii 

12.0 

3.4 
159 ± 38 
104 ± 37 

2.1 

0.3 



TABLE 4: Physical properties of the sample. Halo mass corresponds to the redshift of the arcs; see text for details, 
diameter, AV^, and a listed for J1038 North refer to the entire J1038 system. 
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of the OSIRIS Ha data. We assume a constant stellar 
mass-to-light ratio in the F160W bandpass to calculate 
Sm, ■ We extract the average gas and stellar mass den- 
sities in radial bins along a psuedo-slit oriented along 
the kinematic major axis. (The same slits are later used 
to measure kinematics and metallicity gradients.) We 
examined the gas fractions as a function of radius and 
found it to be nearly constant (< 15% variation). 

2.6. Emission Line Fitting 

The analysis presented in later sections is based on the 
resolved properties of emission lines in the galaxy source 
plane. We determine line flux, velocity, and velocity dis- 
persion from Gaussian fits to the emission lines of in- 
terest. The emission line fit s for J1206 are described in 
detail in I Jones et alj ()2010b[ ). We follow identical meth- 
ods for the remaining sources presented here, except that 
emission lines are fit in the source plane rather than the 
image plane. Here we describe the complete process. 

The OSIRIS data for each object is aligned with high 
resolution HST images used to define the lens model 
transformations from image to source plane f H2.2p . Each 
data cube is reconstructed in the source plane and 
smoothed with a Gaussian kernel to increase signal-to- 
noise of the fainter emission lines and in regions of low 
surface brightness. The smoothing kernels have FWHM 
~ 300 pc for J0744 and J1038, and 600 pc for J1148. We 
then fit a Gaussian profile to the Ha line at each spatial 
pixel using a weighted to account for increased noise 
at the wavelengths of bright sky lines. Weights are deter- 
mined from the noise measured in a blank sky region of 
each data cube. Line fiux and velocity are derived from 
the centroid and area of the best fit profile in all cases 
with signal-to- noise > 5. Intrinsic velocity dispersion 
is calculated by subtracting the instrumental resolution 
ainst — 50 km s~^ (measured from bright OH sky lines) 
in quadrature from the best-fit line width. 

We fit all other emission lines of interest using the ve- 
locity and dispersion derived from Ha as constraints. 
Ha is used because it has the highest signal-to-noise in 
all cases. The centroid and width of each line is fixed 
at the best-fit values measured for Ha, and a weighted 
fit is used to determine the normalization. Each emis- 
sion line is also fit separately without fixing the centroid 
and width, and no significant differences are found. We 
therefore use the emission line fluxes derived with Ha 
line profiles since these have lower formal uncertainties. 
Additionally, when considering the ratio of various emis- 
sion lines, we can be confident that ratios are not biased 



by different kinematic structure in different lines. 

We estimate the physical resolution from observations 
of tip/tilt reference stars taken at the same time as the 
emission line data. Tip/tilt stars are reconstructed and 
smoothed in the same manner and the resulting FWHM 
is measured along the major and minor axes of recon- 
structed star images. These values are reported in Ta- 
ble [TJ Although smoothing degrades the resolution, the 
FWHM is < 1 kpc for all sources. 

3. KINEMATICS 

In order to understand the variation in metallicity gra- 
dients we see in our sample, it is useful to characterize 
each source in terms of its kinematic properties. The de- 
gree of ordered rotation has been determined from fits to 
the Ha emission line ( ij2.6p . The kinemati c properties 
of J074 4 and J1206 are a l ready d escribed in lJones et al.l 
(|2010al) and iJones et"al] (|2010bf) respectively. Here we 
follow a similar analysis for the remaining sources. As 
in our earlier work, we extract one-dimensional velocity 
and dispersion profiles in a pseudo-slit oriented along the 
kinematic major axis (i.e. the direction of highest veloc- 
ity shear). These profiles are shown in Figure [3] along 
with the two-dimensional velocity maps of each source. 
Three galaxies (J0744, J1148, and J1206) show ordered 
rotational motion in both the 2-D maps and 1-D pro- 
files with high local velocity dispersion > 20 kms^^in all 
cases. The fourth, J1038, is comprised of multiple spa- 
tially and kinematically distinct regions indicating that 
this source is undergoing a merger. The diameter and 
peak-to-peak velocity shear AV of each source (measured 
along the pseudo-slit) are given in Table S) For the ro- 
tating sources, AV is related to the maximum circular 
velocity Vmax and inclination angle i as 

AV = 2Vrr,axSilli. (6) 

Additionally we give the mean local velocity dispersion 
(7, defined as the unweighted mean of individual pixels 
with error bars refiecting the la scatter. We note that 
the adopted definition of a varies throughout the relevant 
literature. 

We now briefly compare the kinematics with other 
samples at high redshift reported in the literature. 
The largest such sample is the SINS survey; observa- 
tions of 80 SINS galax i es at z ^ 2 are described in 
iForster Schreiber et al.l (|2009f l. The median velocity 
shear of their sample is AV = 134 kms~^with an in- 
terquartile range of 90 — 220 kms^^. Approximately one 
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third of the SINS sources are clasified as mergers while 
the re st have varying degrees o f rotation and random mo- 
tion. iForster Schreiber et al.l (|2009D tabulate the ratio 
of velocity shear to velocity dispersion, lS.V/{2a), which 
gives an indication of the degree to which each galaxy is 
dynamically supported by rotation vs. random motions. 
The SINS sample has a median AV^/(2(t) — 0.56 and in- 
terquartile range 0.37 — 0.75. Other published samples 
of z ~ 2 — 3 gala xies have r nedian /SV/{2a) va lues of 0.4 
(|Law et al.l 1200 9') and 0.9 (| Jones et al.ll2010a[ ). In com- 
parison, velocity shear observed in the lensed galaxies 
presented here is somewhat higher than in other samples 
at similar redshift. The values in Table 2] give a range 
Ay/ (2it) = 0.76 — 1.4 for our lensed galaxies, and we note 
that our values of a are systematicall y highe r than for the 
definition used by Forstcr Schreiber et al.l (j2009( ). The 
lensed galaxies therefore have a higher degree of rota- 
ti qn than typical SINS sou r ces. A ccording to the criteria 
of IForster Schreiber et all (|2009[ ). all three non-mergmg 
sources in our sample are "rotation-dominated" (defined 
as having AV^/(2ct) > 0.4). 

A common hypothesis advanced by recent studies is 
that the evolution of star-forming galaxies at high red- 
shift is driven i n large part by global gravitational in- 
stability (e . g. I Jones et al.l l2010at iGenzel et all 120081: 
iDekel et al.ll2009D . Briefly, if the velocity dispersion and 
rotational velocity are too small, the system is unstable 
and gravitational perturbations will grow exponentially 
into giant "clumps" . This is quantified by the Toomre 
parameter Q, where values Q < 1 indicate that a galaxy 
is unstable and will fragment into giant clumps. We 
calculate Q/ sini for the rotating galaxies in our sam- 
ple using the inf ormation in Table |4| and following the 
same method as iJones et al.l (|2010a[ ). We list the val- 
ues in Table in as well as the dynamical mass estimated 
as Mdyn — 5R<J^ /G. The results indicate that all three 
galaxies are gravitationally unstable unless the inclina- 
tions are very high (> 60 degrees), althoug h the uncer- 
tainty is approximately a factor of 2 in Q (see Jone s et al.l 
[2010a). This can explain the "clumpy" morphologies ap- 
parent in Figure |2l 

In summary, three galaxies in our lensed sample are 
rotating and one is undergoing a merger. This merger 
fraction is consistent with larger samples at similar red- 
shift. The kinematics of galaxies in our sample are more 
rotation-dominated than others studied at similar red- 
shift, and all rotating galaxies in our sample are likely 
to be gravitationally unstable. The ratio of shear to ve- 
locity dispersion, V/a, is 2x higher than for typical 
galaxies in the extensive SINS survey. 

4. GAS PHASE METALLICITY 

The multiple emission lines detected by OSIRIS enable 
us to determine spatially resolved gas-phase metallicities 
in each galaxy. Robust metallicity measurements rely on 
knowledge of the temperature and ionization state of a 
gas. However, this requires accurate measurement of di- 
agnostic emission lines which are usually too faint to be 
detected at high redshift with current instruments. In- 
stead we use ratios of strong emission lines which can 
be detected in modest integration times and which cor- 
relate with metallicity through locally-calibrated rela- 
tions. These include ratios such as ^^^^^Iro (the 



R23 index) and [Nil] /Ha (the N2 index) for which 
calibrations have been determined using direct elec- 
tron temperature- based measuremen t s or photoioniza- 
tion models (e.g. iPaeel et al.l Il979l: iKew lev & Dooit^ 
'2002: Pcttini fc Pagell l2004HMaiolino et al.ll2008i) . As in 
Jones et al. (2010b) we primarily use the N2 index to 
determine gas-phase oxygen abundance via 

12 + log O/iJ = 8.90 + 0.57 x log [Nil]/ Ha (7) 

(jPettini fc Pagell 1200-5 ). with a dispersion of 0.18 dex. 
The main advantage of N2 is that the two requisite lines 
are very close in wavelength such that systematic un- 
certainties from reddening and instrumental effects are 
negligible. 

While the N2 index provides a practical means of esti- 
mating metallicity in high-redshift galaxies, it can fail in 
cases where {i) active galactic nuclei (AGN) or shock 
excitation contribute significantly to the emission line 
flux, (m) secondary production of nitrogen leads to vari- 
ations in the N/0 ratio, or (m) [Nil] cooling saturates 
(at 12 + logO/H > 9.0). This last case is not a concern 
for the galaxies in our sample as they have metallicities 
significantly below this value. However AGN, shocks, 
and variations in the N/0 ratio are a potential problem 
and we address this issue with additional observations of 
[Oiii] and H/3 emission lines. 

AGN and shocks can be distinguished from star- 
forming regions on the basis of emission line fiux ra- 
tios. One of the most widely used diagnosics is the ra- 
tio of [Nil] /Ha comp ared to [Oiii]/H/3, as described by 
Baldwin et alj (fl98ll ) (the "BPT diagram"). We show 
the BPT diagram for individual pixels in each lensed 
galaxy in Figure j?] Only regions where all requisite 
lines are detected (> Sti) are shown on the diagram. We 
note that some individual pixels are correlated since the 
point spread function is larger than the pixel size. In the 
case of J1206 we estimate H^ flux fro m the global ratio 
Ha/H/3 = 4.07 ()Hainline et al.ll2009l ): these [Oiii]/H^ 
values are therefore accurate on average but may devi- 
ate significantly within individual pixels. The line ratios 
of ~122,000 galaxies from the Sloan Digital Sky Survey 
(SDSS) are also plotted and clearly show a locus of star 
forming galaxies with a separate branch of AGN. The 
lensed galaxies generally follow the locus of local star 
forming galaxies but are offset to higher [Oiii]/H/3 ratios; 
this has been commonly observed in high redshift galax- 
ies and is generally attributed to high ioi iization param- 
eter s rather than AGN (s ee discussion in lErb et al]|2010l 
and iHainline et al.l I2009D . Some regions of the lensed 
galaxies have [Oiii]/H/3 formally above the theoretical 
limit from star formation (Kewley et al. 2001), but con- 
sistent within the uncertainties. We conclude that emis- 
sion lines in the lensed galaxies originate predominantly 
from star forming Hii regions with little contribution 
from AGN or shocks. 

To test whether metallicities derived from the N2 in- 
dex suffer from systematic errors due to variable N/0 
abundance or any other effect, we calculate metallicities 
using various other methods. Using the same data as for 
the BPT diagram (Figure d, we compute 1 2 -flog( 0/H) 
using the cali b ration s of " IPettini fc Pageil (|200l and 
iMaiolino et al.l (|2008[) and show the results in Figure |5j 
All metallicities are consistent given the intrinsic scatter 
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FIG. 3: Source plane kinematics of the lensed galaxies. Left: two-dimensional velocity field. Right: one-dimensional velocity 
(red circles) and dispersion (black triangles) of each source, extracted along the kinematic major axis. 
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FIG. 4: Diagnostic diagram of [Nn]/Ha and [Olll]/H/3. 
Eacli large point represents a single spatial pixel of OSIRIS, 
and the small grey points are ^^122, 000 galaxies from the 
Sloan Digital Sky Survey (SDSS) with signal to noise > 5 in 
all relevant emission lines. SDSS galaxies show a locus of star 
forming galaxies, and a separate branch of AGN extending 
to the upper right. The dotted line shows t he theoretical 
maximum [Olll]/H/3 from star forming regions (jKewlev et al.l 
I2001D . and the dashed line is an e mpirical division betwee n 
star-forming galaxies and AGN from lKauffmann et al.l ()2003l ). 
All regions of the lensed galaxies are consistent with pure star 
formation, although [Olll]/H/3 is typically above the locus 
of SDSS star forming galaxies at z ~ 0. This is commonly 
observed at high rcdshift and is attributed to a high ionization 
parameter (e.g. Mainline ct al. 2009; Erb ot al. 2010). 

in each calibration, although [Oiii]-based metallicities are 
typically lower by ^ 0.15 dex than those inf e rred f rom 
[Nil] alone. Furthermore the IMaiolino et al.l (I2008D re- 
sults are systematically higher than for IPettini fc Pagell 
( j2004i) by ~ 0.2 dex. Overall we conclude that systematic 
uncertainties in metallicity calculated from the N2 index 
are limited to < 0.2 dex depending on the calibration 
used (e.g. N2 vs. 03HB), with an addi tional ^0.2 dex 
uncer taint y in the zero point (e .g. from IPettini fc Pagell 
120041 vs. iMaiolino et"all 120081 ). In the following sec- 
tion we show how different calibrations affect the derived 
metallicity gradients. 

4.1. Metallicity Gradients 

We are now in a position to determine gradients in the 
inferred metallicity along the kinematic major axis of 
each galaxy. To do so we extract resolved measurements 
of [Nil] and Ha flux along a pseudo-slit in the same 
manner as for constructing the rotation curves shown in 
Figure [3] While other orientations can also be used, they 
must be corrected by the uncertain inclination and so we 
take the major axis to provide the most precise measure- 
ment. In the case of J1038 the inclination angle between 
merging components is unknown, hence the magnitude of 
gradients for this source should be considered an upper 



limit. The center of each galaxy is defined kinematically 
as the position along the pseudo-slit at which the veloc- 
ity is equal to the systemic velocity (i.e. velocity = 
in Figure EI; for the merging system J1038 we use the 
systemic velocity of the UV-bright northern region (cor- 
responding to velocity ~ 100 kms~^ in Figure |3]). We 
bin individual pixels by their galactocentric radius and 
compute the total Ha and [Nil] flux using all pixels with 
signal-to-noise > 10 in Ha. Metallicity is derived from 
the resulting [Nil] /Ha ratio via Equation [7] and shown 
as a function of radius in Figure |6l 

Each galaxy shows significant variations in metallicity 
as a function of radius (Figure [51). The three rotating 
galaxies (J0744, J1148, and J1206) all show significant 
gradients with decreasing metallicity at larger radii, as 
observed in all local disk galaxies. The merging system 
J1038 has significantly lower nuclear metallicity than the 
rotating galaxies, and exhibits a more complex metallic- 
ity distribution. Metallicity is measured for three spa- 
tially distinct regions of J1038. The peak of rest-UV 
and Ha emission (at i? = in Figure [5]) has the low- 
est metallicity at 12 + \ogO/H = 8.2, while the peak of 
IRAC 3.6/im flux (at i? = 5 kpc) has marginally higher 
12 -I- \ogO/H = 8.4. The third component, located be- 
tween the first two (at i? = 3 kpc) has the highest mea- 
sured metallicity. 

We now quantify the radial metallicity gradient of each 
source in order to compare the results with other sam- 
ples. Gradients in local galaxies are commonly expressed 
as a linear relation in units of dex/kpc. We compute 
the best linear fit to 12-1- \ogO/H as a function of ra- 
dius using a weighted least-squares method. The best 
fit relations are shown as dashed lines in Figure [5] (right 
panels) and provide a reasonable fit to the data with re- 
duced values of 0.3 — 2.3, although we note that the 
linear fits overpredict metallicity at the largest radii in 
all sources. The best-fit central metallicity and gradient 
of each source is listed in Table [S] 

To examine whether these results might be affected 
by systematic calibration errors, we compute the metal- 
licity gradient in the exact same manner using different 
strong-line metallicity diagnostics. In Figure [7] we show 
the radial metallicity profiles determined from all a vail- 
able calibrations described bv lPettini fc Pagell (|2004j ) and 
Maiolino ct al. (2008) as well as their best-fit linear gra- 
dients. We note that J1206 is assumed to have a con- 
stant Ha/H/3 ratio, that [Oiii] and H/3 in J1038 were 
only observed in the region corresponding to i? < 4 kpc, 
and that [Oiii] and H/3 are only detected with sufficient 
signal-to- noise at i? < 1 kpc in J1148. We further assume 
a constant Ha/H/3 ratio for J0744 since H/3 is detected in 
only a few pixels. Systematic offsets between the various 
calibrations are apparent in Figure [71 but in general the 
different methods agree given their intrinsic scatter. We 
give the best-fit gradients in Table [3 In most cases the 
resulting gradients are stronger than that for our adopted 
N2 calibration, and in all cases the uncertainty is larger. 
Importantly, [Oiii]-based calibrations (including 03HB 
which is independent of N2) confirm that line ratio gra- 
dients are due to metallicity gradients rather than N/0 
abundance or other effects. The exception is J0744 for 
which [Oiii] and H/3 ratios suggest a positive metallic- 
ity gradient. Overall, the different metallicity derivations 
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FIG. 5: Comparison of the different strong-line metallicity diagn ostics N 2 ([Nllj/Hia), 03N2 ([Oni]/[ Nll]), and 03HB 
([Oin]/H^). Metallicities are derived from the calibrations of Pettini fc Page] (j2004) (left panel) and Maio lino et al.l (|2008l ') 
(center and right p anels). Each point rep resents an individual OSIRIS pixel; data are identical to those in Figure 2] The la 
scatter measured bv lPettini fc Page] (|2004f ) is shown in the lower right of the left hand panel. Metallicities derived from different 
methods are generally i n agreement, w i th [Ol llj-based metallicities typically 0.15 dex lower than those derived from N2. Metal- 
licity derived using the lMaiolino et al] l|2008l ) calibrations are systematically higher by ~ 0.2 dex compared to lPettini fc Pagell 
l|20M ). 



are in reasonable agreement, albeit with large uncertain- 
ties due to both low signal-to-noise and intrinsic scatter 
in the calibrations. In the following sections we will con- 
sider the N2-based gradients of each galaxy to be the 
most accurate. 

4.2. Evolution with Redshift 

The metallicity gradients we have determined for our 
lensed sample can now be compared to those for galax- 
ies at different redshifts. Figure [8] shows a comparison 
with other measurements reported in the literature. It is 
instructive to separate the sample according to whether 
the galaxies are isolated or interacting. Figure [5] demon- 
strates that interacting galaxies tend to have flatter gra- 
dients and lower central metallicity. Two of the lensed 
z 2 sources, J1148 and J1206, have gradients < —0.2 
dex/kpc which are significantly steeper than those found 
in local massive disk galaxies. The only local galaxies 
in Figure [5] with similar or steeper gradients are barred 
spiral starbu rsts from [Co nsidcrc ct al. (2000) and dwarf 
galaxies from ! Vila- Costas fc Edmundsi (11992) . J0744 has 
a gradient typical of local disk galaxies. The merging 
galaxy J1038 has an "inverted" gradient (> dex/kpc) 
- a phenom enon which has also been observed in lo- 
cal mergers ( Rich et al.ll2012D and high -redshift sources 
(jCresci et al.l 120 lOtlQuevrel et al.ll2012D . 



studied the stellar to halo mass ratio and its variation 
with redshift by combining cosmologic al simulations with 
measured stellar mass fu n ctions (e.g. lYang et al.l 120111 : 
iConrov fc W echsler' '2009; Moste r et al.l I2012D . We use 
the results of Moster ct al. (2012!) to do this and report 
the values in Table [H We calculate the corresponding 
halo mass at different redshifts b y integrating th e mean 
halo growth rate determined by iFakhourT "etaLl dlOlO), 
and determine the stellar mass at these redshifts via the 
iMoster et al.l ()2012D formalism. Next we must determine 
galaxy stellar masses for the comparison samples. For 
the MASSIV sample, we use the stellar masses derived 



Although the overall sample size is still limited, it ap- 
pears that the metallicity gradients vary systematically 
with redshift and here we seek to quantify this evolution. 
However, before doing so, it is important to consider how 
best to compare galaxies at different redshifts in a con- 
sistent manner. One way to proceed is to determine the 
approximate stellar mass that the lensed galaxies will 
have at different cosmic epochs based on halo abundance 
matching techniques. This enables us to select a com- 
parison sample with stellar masses approximately equal 
to the value inferred for the lensed galaxies at that red- 
shift. The first step in this process is to estimate the 
halo mass for our lensed galaxies. Several groups have 



bv lQuevrel et all (|2012D . For aU other comparison sam- 
ples at z ~ 0, we extract broadband B — V color and ab- 
solute luminosity {Mb and/or Mk) from the NASA Ex- 
tragalactic Database and calculate the stellar mass using 
the m ass-to-light ratio formulae from Torres-Flores et all 
(|2011h . Finally, we construct two separate comparison 
samples appropriate for (1) isolated and (2) interact- 
ing galaxies. For isolated galaxies we choose a stellar 
mass range corresponding to = lO^ '^^^-^ a.t z = 2.2; 
the isolated lensed galaxies are all within 0.3 dex of this 
value. The stellar mass range for interacting galaxies 
corresponds to that of J 1038. We show the metallicity 
gradients of each comparison sample as a function of lin- 
ear cosmic time and equivalent redshift in Figure [S] 

Figure [5] shows a clear trend in metallicity gradients 
with time (or redshift) for our sample of galaxies which 
occupy similar dark matter halos. The scatter at a given 
redshift presumably reflects the degree of intrinsic vari- 
ation within the population. To investigate this with 
more clarity, we show in Figure [TO] the mean and la 
scatter of each isolated galaxy comparison sample. Con- 
sidering only the z > 2 and z data in Figure [TOl 
we can see that the average metallicity gradient becomes 
flatter with time although the intrinsic scatter is simi- 
lar to the average gradient. However, data at z ^ 1.2 
from the MASSIV survey do not support this picture; 
instead they show a mean gradient ~ with relatively 
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FIG. 6: Metallicity gradients of the lensed galaxies. Left: [Nil] /Ha ratio as a function of radius. Right: Gas-phase metallicity 
as a function of radius derived from the N2 index calibration of iPettini fc Paeeil (|2004l ). Dashed lines show the best fit linear 
metallicity gradient. [NiiJ/Hq ratios are extracted along the kinematic major axis using the same slit as for the right panels of 
Figure [3l 
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FIG. 7: Metallicity gradients for each galaxy in our sample based on the use of various strong-line d iagnostics. N2, 03N2 , 
and 03HB are as presented in Figure (5] and we use the calibrations of lPettini fc Pagell (gOOl) (PP04) and lMaiolino et al.l (|2008l ') 
(M+08). Solid lines show the gradient from Figure [6] and dotted lines are the best fits to the other diagnostics. For J0744 and 
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Name Central 12 + log O/H N2 PP04 N2 M+08 03N2 PP04 03N2 M+08 03HB M+08 

(from N2, PP04) (dex/kpc) (dex/kpc) (dcx/kpc) (dex/kpc) (dex/kpc) 

MACS J0744 8.67 + 0.18 -0.06 + 0.04 -0.13 + 0.05 0.02 + 0.04 0.02 + 0.04 0.10 + 0.07 

SDSS J1038 8.26 + 0.19 0.08 + 0.03 0.15 + 0.07 0.25 + 0.07 0.31 + 0.08 0.37 + 0.09 

SDSS J1148 8.72 + 0.18 -0.28 + 0.05 -0.51 + 0.11 -0.33 + 0.12 -0.37 + 0.11 -0.37 + 0.19 

SDSS J1206 8.55 + 0.18 -0.25 + 0.06 -0.45 + 0.08 -0.40 + 0.08 -0.46 + 0.08 -0.40 + 0.11 



TABLE 5: Metallicity gradients derived from different strong-line diagnostics shown in Figure [T] 
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low scatter of ~ 0.05 dex. This discrepancy is puzzling 
and warrants further investigation. Formally the mean 
metallicity gradient for the samples shown in Figure [10] 
is -0.20 ± 0.07 at z = 2.2, 0.005 ± 0.011 at z = 1.2, 
and —0.077 ± 0.005 at z = 0. The local data and lensed 
galaxies suggest that on average, metallicity gradients 
have declined by a factor of 2.6 ± 0.9 in the past 10 Gyr. 
Naturally metallicity gradients as steep as —0.2 dex/kpc 
cannot be common in massive galaxies at z ~ simply 
because of their size. In the Milky Way, solar metallicity 
at i? = 8 kpc would imply 40 x solar metallicity in the 
center with such a gradient. 

We now consider the effect of "inside-out" size growth. 
Supposing that the relative metal enrichment at the cen- 
ter and effective radius Re of a galaxy remains constant 
with time, the gradient will scale as 1/Rp, and will flat- 
ten wi th time as the radius increases. Ivan Dokkum et all 
(|2010( ) have measured the size of massive galaxies se- 
lected to have a fixed number density at different red- 
shifts (equivalently, approximately the same halos), and 
find i?e oc (1 -f z)~^'^^. We show the appropriate change 
in metallicity gradient (cx (1 -I- z)^-^^) in Figure [101 and 
find that it is within 1 cr of the Milky Way and z = 
data. However, the galaxies used to determine this size 
growth are significantly more massive (by 0.7 dex at 
z = 0) than the metallicity gradient sample and likely 
grow at a different rate. Adopting the same functional 
form Re oc (l-|-z)~", the lensed galaxi es must grow in size 
with a — 0.8 ± 0.4 (cf. a = 1.27 from Ivan Dokkum et al.l 
|2010( ) to match the z = comparison sample mean. 
Metallicity gradients should therefore flatten with time 
as a natural consequence of inside-out growth, and the 
amount of size growth required is in rough agreement 
with direct measurements. 

5. CHEMICAL EVOLUTION MODEL 

As we have seen a first clear trend of flattening in the 
metallicity gradients for comparable systems with cosmic 
time, we seek now to develop a physically-motivated ex- 
planation of the trend. Our purpose in this section is to 
explain the metallicity of a star-forming system (specif- 
ically a galaxy or a region within a galaxy) in terms of 
gas accretion, outflows, and star formation. The sim- 
plest such model is the "closed box" in which accretion 
and outflow rates are both zero. In this case the gas- 
phase metal mass fraction Z is determined entirely by 
the yield y and gas fraction /gas as 

Z = yln /" s 



(8) 



(|Schmidtl[T96l . This can be equivalently expressed in 
terms of the oxygen number abundance. 



l2 + logO/H = 12 + log 



yo 



11.728 



-flog(ln/-i) (9) 



where yo is the oxygen yield and the factor 11.728 
converts oxygen mass abundance to number abundance 
(|Lee et all [2001 . Throughout this paper we assume a 
value yo = 0.0087, appropriate for solar abundance ra- 
tios with y = 0.02 ijFinlator fc Dave"2008'). 

We now construct a model which includes accretion 
and outflow of gas. The star-forming system which we 
model has the following properties: a gas mass Mg, 
stellar mass M*, star formation rate SFR= M,, total 



mass of metals in the gas Mz, and metallicity defined as 
Zg — Mz/Mg. We consider four factors which affect the 
metallicity at a given time: 

1. Inflowing gas. The gas inflow rate is assumed to be 
a constant multiple fi of the star formation rate. 
Mi = fiM*. Inflowing gas is assumed to be metal- 
free [Z = 0) corresponding to pristine gas accreted 
from the intergalactic medium. 

2. Transported gas. In addition to inflow Mi, we as- 
sume that some accreted gas is already enriched 
with a metallicity Z = Zg and accreted at a con- 
stant multiple of the star formation rate. Men = 
fenM^,. If the system considered is a region within 
a galaxy, then this term corresponds to the trans- 
portation of metal-enriched gas within the galaxy. 
If the system is a galaxy, this can be thought of 
as a "galactic fountain" term which describes the 
re-accretion of gas which has been expunged in out- 
flows. 

3. Outflowing gas. Outflows are a vital component 
in any reasonable model of galaxy chemical evolu- 
tion. Outflows are observed ubiquitously in galax- 
ies with star formation rates > 0.1 MQyr^^kpc^^ 
( HeckmanI 120021 ) which is greatly exceeded by the 
lensed galaxies studied in this paper. The outflow- 
ing gas is assumed to have metallicity Z — Zg and 
a rate proportional to the SFR, Mo — foM^,. The 
ratio of outflow rate to SFR, /o, is called the mass 
loading factor (MLF). 

4. Star formation. Star formation is ultimately the 
source of all heavy elements. The amount of met- 
als produced per unit time is deflned as the yield 
multiplied by the star formation rate, yM^. Star 
formation additionally removes metals from the gas 
reservoir at a rate ZM^. 

We additionally assume that the gas mass remains con- 
stant. With these definitions, the metallicity is expressed 
as 

2/0 



12 + log O/i? =12 + log 
+ log 



11.728 
l-exp( -(1 + /^) 



log(l + /^) 



/ga 



(10) 



where we define the "effective mass loading factor" 
(EMLF) /o = /o — fen- (A complete derivation of Equa- 
tion [10] is provided in Appendix [X] along with an assess- 
ment of the assumptions made.) Therefore, our model 
provides the gas phase metallicity as a function of EMLF 
and gas fraction. Equivalently, can be determined 
from measurements of metallicity and gas fraction. We 
show the results of this model for various values of in 
Figure [TT] along with measurements of the lensed galax- 
ies. 

5.1. fo as the Origin of Metallicity Gradients 

The presence of metallicity gradients in local disk 
galaxies is typically attributed to variations in the mass 
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FIG. 11: Gas phase metallicity as a function of gas fraction. 
We show predictions from our chemical evolution model as 
well as the closed box model. Data points correspond to 
measured values of the galaxies J0744, J1148, and J1206. 

loading factor with radius and/or inward radial migra- 
tion of metal-enriched gas. Both of these physical effects 
can be expressed as a variation of with radius, where 
negative metallicity gradients imply that increases 
with radius. In Appendix |^ we derive the relation be- 
tween gradients in metallicity and f^, and here we apply 
these findings to the lensed galaxies. Equation lAlll shows 

that a gradient in (expressed as Ar = ^) naturally 
gives rise to a metallicity gradient Az = -^(\ogO / H). 
For fiducial lensed galaxy values /gas — 0.5 and ~ 2, 



-{).12Ar 



dex. A gradient Ai^ = 0.5 — 2 kpc^^ can therefore explain 
the observed metallicity gradients of the isolated galaxies 
(Table [5]), and the merger J1038 requires An = —0.7 
kpc~^ to explain the inverted gradient. 

Let us now consider the evolution of such a gradient 
with time under the assumption that the gas fraction 
does not vary significantly with radius. In this case the 
time derivative is given by Equation IA12I An im- 
portant qualitative feature of this equation is that 

is negative if > 0, in which case metallicity gra- 

dients will become stronger with time. This is easily 
understood in terms of Figure 1111 Physically, we have 
assumed that the gas fraction is constant throughout a 
galaxy and declines monotonically with time. As fgas de- 
creases, the difference in metallicity between two regions 
with different will become larger and hence gradients 
will become more pronounced. However, it is unlikely 
that An increases with time. For example, taking a 
typical value = 2 at a radius R = 1 kpc, the even- 
tual metallicity is only 12 + \ogO/H = 8.4 or about 
half the solar value. Descendents of these lensed galax- 
ies typically have super-solar metallicity at such small 
radii, implying that must decrease at later times. The 
required decrease in mass loading factor can be quanti- 
fied from the observational constraint that < 0.05 
dexkpc"^ Gyr~^ (Figure ITU)) . Using fiducial values for 
the lensed galaxies /gas — 0.5, /^ ~ 2, Af* ~ 100 
Mgyr-i, Mgas =i 2 x 10^° Mq, and Ai^ ~ 2 kpc-\ Equa- 
tion Ell] implies that ^ ~ -4 kpc-iQyr-i. There- 
fore the effective mass loading factor must decrease on 
a time scale ti = A^ (^^)~^ 500 Myr, compara- 
ble to the characteristic age of recent star formation in 
the lensed galaxies. This result suggests that a phase of 
strong starburst-driven outfiows with high mass loading 
factors lasts only a few hundred Myr, which matches the 
duty cycle of strong star for mation in LBG s inferred from 
statistical studies (e.g. Stark et al.ll2009D . 

In summary, we have used the chemical evolution 
model to show that a radial gradient in the effective mass 
loading factor Aji produces a metallicity gradient. The 
observed time evolution of metallicity gradients requires 
that An must decrease with time, and the characteris- 
tic starburst-driven outfiow time scale from our model 
(~ 500 Myr) is in rough agreement with the duty cycle 
of star formation inferred from statistical measurements 
of high redshift LEG demographics. 

6. DISCUSSION 

We now discuss what can be learned physically about 
the mass assembly in z ~ 2 galaxies by applying the 
model discussed in fj5]to measurements of the galaxies in 
our sample. 

6.1. Mass Loading Factor 

The mass loading factor, defined as the ratio fo of mass 
outflow rate to star formation rate, is an important ingre- 
dient of galaxy formation models. Our simple chemical 
evolution model can be used to infer fo from measure- 
ments of metallicity and gas fraction via Equation 1101 
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We apply this equation to the integrated values of the 
lensed galaxies (Table |4]) and list the resulting effective 
mass loading factors in Table [6] Assuming that ac- 
creted gas has an average metallicity significantly lower 
than that of the galaxies, we s et /en ~ and equate 
f'o — fo (see discussion in i 3A.l| . Therefore, the model 
implies a modest mass loading factor fo < 2.3 for all 
three galaxies studied in detail here. 

In the case of J1148, we can estimate the true mass 
loading factor from direct measurements of the outfiow 
column density. Quider et al . (2009) measure a column 
density logiVsin = 16 cm~^ in the outfiowing gas. The 
Siii/Siiv ratio is ~ 12 suggesting that most of the out- 
flowing mass is in neutral Hydrogen gas traced by Sill. 
Motivated by our chemical evolution model and the re- 
sults discussed later in this section (cf. Figure [T^ , we 
assume that most of the outfiow occurs at large ra- 
dius with metallicity 12 -|- logO/H ~ 8.1 (Figure |6]) or 
[0/H] ~ —0.6. Since O and Si are both a elements gen- 
erated by the same nucleosynthetic processes, we assume 

that [Si/H] 0.6 as weU. Indeed, iPettini et all ()2002l ) 

show that variation in a element abundances is < 0.2 dex 
in the ISM of the Lyman break galaxy cB 58. Adopting 
a sol ar abundance ratio [O/H]© = —4.44 (jPettini et al.l 
[200l gives an outfiow column density log Nh = 21.0 
cm~^. This corresponds to a total mass column density 
(including Helium) = 1.1 x 10'' MQkpc"^. We can 
crudely estimate the total mass assuming that outflowing 
gas uniformly fills a sphere of radius R: 

Mo = ^TrR^T.o - 3 X 10^ Mq 

for R = AwM^/M* ~ 8 kpc, outflow velocity Aw = 200 

kms-^ ()Quider et al.l [20091 ) . and timescale M^/M^ = 40 
Myr (Table ini). The time-average mass loading factor is 
then given by fo = Mo/M^, ~ 0.4. Uncertainty in 
due to the assumed geometry is approximately an order 
of magnitude, such that we can confidently assert that 
/o < 5 in reasonable agreement with the value derived 
from our chemical evolution model. 

We now turn to spatial variations in the mass loading 
factor, and radial gradients in particular. Gas fractions 
in the lensed galaxies have almost no variation with ra- 
dius, and therefore we can apply the formalism devel- 
oped in i jS.ll to estimate /S.r, and the characteristic 
timescale ti. The results are listed in Table [6l We also 
give the mass assembly timescale Mtf/M^,. We can com- 
pare the simple analytic treatment of §5.1l with variations 
in f'o inferred directly from spatially-resolved measure- 
ments via Equation [TU] We use the measured radial pro- 
files of gas fraction and metallicity to compute /q, and 
show the results in Figure [121 As expected, f'o increases 
at large radii. The analytic approximation is in rea- 
sonable agreement with a linear fit to the data, although 
the measured radial gradient of f'o is poorly represented 
by a straight line. Metallicity gradients observed in the 
lensed galaxies are therefore likely to originate from ra- 
dial variations of the effective mass loading factor. 

6.2. Inflow Rate 

Under the equilibrium condition that the gas mass re- 
mains constant, the infiow rate is directly related to star 
formation rate and mass loading factor via Equation I A31 



When considering a galaxy as a whole, /en is negligible 
and thus /,; = 1 + /o, or Mi = (1 + fo)M^. We calculate 
Mi using the relevant values in Tables [4| and [6| and list 
the results in Table [51 Infall rates range from 1.2 — 2. Ox 
the star formation rate. 

6.3. Radial Gas Transport 

We have suggested that radial metallicity gradients are 
caused by gradients in the effective mass loading fac- 
tor. Negative gradients can be caused by increasing mass 
loading factor with radius (^ > 0), by radial transport 

of enriched gas (^f^ < 0), or both. The degeneracy be- 
tween fo and fen captured by Equation IA5I essentially 
means that the chemical evolution model cannot dis- 
tinguish the contributions of outfiow and gas transport. 
However, in cases where /^ < 0, gas transport is neces- 
sarily present. Figure [T^l shows that the inner regions of 
J 1148 are inferred to have /o < with 2cr significance 
from uncertainty in the gas fraction, although the signif- 
icance is marginal when considering a number of other 
assumptions made (e.g. the Schmidt relation, constant 
mass-to- light ratio, and the model itself). Nonetheless 
the high metallicity and high inferred gas fraction suggest 
that gas in the inner regions of J 1148 has likely been en- 
riched by previous star formation episodes at larger radii. 
This observation, combined with the clumpy structure 
and low Q < 1 of J 1148, are in qualitative agreement 
with the model of iDekel et al.l (|2009| ) in which gravita- 
tional instability leads to the formation of giant clumps 
which migrate into the galactic center in < 500 Myr. 

We reiterate that this result depends crucially on 
the gas fraction, which we have estimated using the 
Kennicutt-Schmidt law. Future measurements of the gas 
density e.g. with ALMA are needed for direct verifica- 
tion. 

6.4. Positive Metallicity Gradients 

In the local universe, disk galaxies almost invariably 
show negative metallicity gradients. Only a few posi- 
tive gradients have been measured, all of which are rel- 
atively shallow (< 0.05 dex/kpc) and found exclusively 
in merging systems (see Figure [5]) • In contrast, positive 
metallicity gradients have been claimed for some high 
redshift galaxies, including many with no sign s of recent 
interaction or merging (Cresci ct al. 2010; Qu evrel et al.l 
l2012f ). In the lensed sample we find that all three iso- 
lated galaxies have negative gradients, while the merging 
system J1038 has a positive gradient. What could ex- 
plain these claimed positive gradients at high redshifts? 
Cresci et al. (2010) and Qucyrcl et al. (2012) have sug- 
gested they could be caused by high accretion rates of 
metal-poor gas which occur only at high redshift or dur- 
ing mergers. Even so it is difficult to form a positive 
gradient since the fraction of gas converted into stars 
per unit time must increase with radius, while the dy- 
namical time scale generally decreases with radius. The 
corollary is that a higher fraction of metal-enriched gas 
must be lost to outfiows at small radii, or in terms of our 
chemical evolution model, fo must decrease with radius. 
While gaseous infall is notoriously difficult to detect at 
high redshift, outfiows are relatively easy to detect and 
in principle fo can be estimated as a function of radius 
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following the method of ij6.1l Current planned and on- 
going observations will directly address the variation of 
fo with radius in high redshift galaxies and determine 
whether this is a viable mechanism of producing positive 
mctallicity gradients (see ij6.5p . 

Alternatively, it is possible that many of the positive 
(and negative) gradients reported at high redshift are a 
result of systematic errors rather than actual metallicity 
gradients. In this paper we have primarily used the ratio 
of [Nil] /Ha as a proxy for metallicity, however this alone 
is not convincing because of the issues discussed in 21 
For example, IWright et alj (|25i0i ) describe a galaxy with 
strong radial line ratio gradients which mimic a negative 
metallicity gradient in [Nil] /Ha and a positive gradi- 
ent in [Oiii]/H/ 3. Using the diagno stic diagram of Fig- 
ure HI however. IWright et al.l ()2010( ) show that these are 
caused by a central AGN surrounded by an extended 
star- forming disk rather than a metallicity gr adient. In 
the local universe, iWestmoauette et"al] ()2012|) find that 
[Nil] /Ha increases with radius in 7/18 galaxies in their 
volume-limited sample and show, via measurements of 
[Sii] and [Oi] emission, that these are caused by an ion- 
ization parameter gradient. The most likely explanation 
is an increasing contribution of shocks at large radius 
rather than a positive metallicity gradient. Therefore 
it is known that the line ratio signatures of both posi- 
tive and negative metallicity gradients are mimicked by 
AGN and shocks, and we caution that gradients based on 
strong line ratios are not reliable unless confirmed with 
multiple appropriate diagnostics (e.g. Figures |4] and [7|) . 

6.5. Future Work 

The study of galaxy formation through time evolution 
of metallicity gradients is still in its infancy. Here we 
briefiy discuss three areas of future work that will im- 
prove our understanding: the mass loading factor, sam- 
ple sizes and redshift range, and gas fractions. 

We have shown from simple chemical evolution argu- 
ments that metallicity gradients likely originate from ra- 
dial gradients in the mass loading factor fo- Figure [12] 
suggests that this quantity is several times higher at large 
radius (i? > 1 kpc for the lensed galaxies) than at the 
center. Furthermore, we have shown that a radially de- 
creasing fo can create the positive metallicity gradients 
observed at high redshift but for which there are no local 
analogs. While difficult to measure accurately, can be 
estimated from rcst-UV absorption lines as described in 
§6.11 and spatially resolved measurements can constrain 
its variation with radius. We are actively pursuing these 
measurements in lensed galaxies at z ^ 2 via an approved 
program with the Keck II/ESI integral field unit. 

The sample of high redshift galaxies with robust metal- 
licity gradient measurements remains very small. Al- 
though data exist for more than 30 galaxies at z > 1, 
the majority lack the data needed to distinguish true 
metallicity gradients from the signatures of shocks and 
AGN. Furthermore, different samples reported at z > 1 
(|Cresci et al.ll20Tol : lYuan et al.ll201lt iQuevrel et aIll20Tl 
this work) are discrepant and warrant further investiga- 
tion. Additionally the range of redshifts 0.1 < z < 1 
has not yet been explored, and our data suggest that 
gradients should decrease in magnitude by a factor of 
~ 2 in this time (e.g. Figure [T0|. It is important to re- 
alize that highly multiplexed optical spectrographs can 



efficiently measure resolved emission line fluxes to large 
galactocentric radius in s tatistically significant s amples 
of galaxies out to z < 1.7 (jMiller et al.ll2011ll2012[ ). This 
is an attractive means of procuring a large number of 
metallicity gradient measurements at intermediate red- 
shifts z < 1, especially in tandem with multi-object near 
infrared spectrographs to detect Ha and [Nil] at z > 0.5, 
and we will exploit these capabilities in future work. 

Finally, results from chemical evolution modelling are 
sensitive to the gas fraction which is difficult to measure 
at high redshift. In the present work we have estimated 
gas masses by inverting the Kennicutt-Schmidt relation, 
introducing a considerable degree of uncertainty. Gas 
mass is difficult to determine even with direct measure- 
ments and typically relies on an uncertain factor Xco to 
convert a CO line luminosity to gas mass. With the ad- 
vent of ALMA it is now possible to resolve the molecular 
gas emission in typical star forming galaxies at high red- 
shift, and calibrate gas mass-to-luminosity ratios from 
independent dyna mical measurements (e.g. iDaddi et al.l 
I2OI0I : iStark et al.l[2008i) . With such advances, we may 
be able to confidently determine resolved gas fractions in 
the near future. 

7. SUMMARY 

This paper presents robust spatially resolved metallic- 
ity measurements of four lensed galaxies at z = 2 — 2.4 
based on the strong emission lines Ha, [Nil], [Oiii], and 
H/3. The combination of gravitational lensing and adap- 
tive optics provides a source plane resolution of up to 
300 — 600 parsecs in each galaxy. Three targets are iso- 
lated rotating galaxies, and one is a merger of at least 
two systems with mass ratio (6 ± 3) : 1. All three ro- 
tating galaxies have negative radial metallicity gradients 
inferred from the [Nil] /Ha emission line ratio and con- 
firmed with [Oiii] and H/3 measurements. The merging 
system has a positive gradient in the UV-bright source, 
also confirmed with multiple emission line ratios. 

Metallicity gradients measured at z > 2 are compared 
with galaxies at lower redshift selected to occupy equiv- 
alent dark matter halos. On average, gradients in the 
rotating galaxies must flatten by a factor of 2.6 ± 0.9 
between z = 2.2 and z = 0. This factor is in rough 
agreement with size evolution mea s ured for more mas- 
sive galaxies by Ivan Dokkum et all (|2010D , hence radial 
inside-out growth can account for the variation in metal- 
licity gradients with redshift. 

We develop a simple model of chemical evolution to ex- 
plain the emerging body of data. From a combination of 
spatially resolved gas fraction and metallicity, the model 
implies that metallicity gradients are caused by radial 
variation in the effective mass loading factor. The time- 
averaged mass loading factor is inferred to be < 2. Fi- 
nally, inward radial gas transport is required to explain 
the high metallicity and gas fraction of J 1148, and is 
likely present in the other galaxies although we cannot 
determine this directly from the model. 

The study of metallicity gradients at different redshifts 
is still relatively new and robust sample sizes are lim- 
ited. Future work addressing direct measurements of the 
molecular gas content, metallicity gradients at intermedi- 
ate redshifts (0.2 < z < 1), and direct, spatially resolved 
measurements of the mass loading factor are planned or 
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TABLE 6: Properties inferred from the chemical evolution model. Effective mass outflow rates and uncertainties correspond 
to the formal uncertainties in gas fraction for the tabulated value of 12 -I- \ogO/H. Ar, its derivative, and ti are calculated 
from /o except in the case of J1148 where we use fo + o' = 1.2. 
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FIG. 12: Effective mass loading factor vs. radius. is calculated from Equation [10] and is shown to increase with radius. 
The measured gas fraction and metallicity of J1148 is unphysical at _R < 1 kpc in the context of our chemical evolution model, 
and so we also show the results of perturbing the gas fraction by multiples of the uncertainty o-(/gas) — 0.15. We note that the 
overall shape and normalization of fo{R) is relatively constant over a range A/gas = 4a = 0.6; clearly Equation [TU] is relatively 
insensitive to gas fraction in this regime. Dotted lines show the result from the simple analytic approximation described in ^5.1\ 
(see discussion in the text). For J 1148 the analytic prediction applies to the case /gas — o". 



ongoing. These new observations will confirm (or refute) 
the results of this paper and will place additional con- 
straints on the formation of galaxies from the origin and 
evolution of metallicity gradients. 
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APPENDIX 

A. CHEMICAL EVOLUTION MODEL 

In this Appendix we provide the mathematical details of the chemical evolution model discussed in the text. All 
notation and definitions are the same as in iJS] To begin, the combination of infiowing gas, transported gas, outfiowing 
gas, and star formation as described in iJ5] gives a rate of change of metal mass 

= yAU - ZgNU - ZgMo + ZgMen. (Al) 

The gas mass is assumed to be constant, i.e. 

Mg = M, + Men ~Mo-M*=0 (A2) 

or equivalently 

h + fen-fo-1^ 0. (A3) 
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Equation I All therefore simplifies to 



yAU- ZgM,~ Zgf^M,. 



where we define the "effective mass loading factor" (EMLF) as 

fo — fo fen • 

The rate of change of metallicity is 



dZg 

dt 



d Mz 



1 dM,_ 
Wg dt 



from Equation IA4I Integrating with respect to time gives 

^ jln{y-Zg-Zg,0 + 



1 + /^ 



1 + /^ 



7ln(y) 



M„ 



gas 



(A4) 
(A5) 

(A6) 
(A7) 



gas 



where we introduce the gas fraction defined in Equation [5] and apply the boundary conditions Zg ~ and Af* = at 
time t = 0. With some algebra we obtain the desired expression for metallicity, 



7 — " 



l-CXp(-(l + /^)-^ 



(A8) 



While Equation IA8I contains all of the information needed, it is more convenient to express the metallicity in terms of 
12+log(0/H). The equation then becomes 



12 + log C'/i7 = 12 + log 



yo 



11.728 



log (1 + /^)+ log 



exp 



/gas 



(A9) 



as given in the text. 

We now derive the relation between gradients in metallicity and and their evolution with time. From Eguation lAQi 
metallicity varies with as 

— log OH) = '\\ AlO 

df'o ln(10)(l + /^)[exp(a;) - 1] 

where a; = (1 + /o)(l — /gas)//gas- This gives rise to a metalhcity gradient of magnitude 



Az = -^{logOlH) = An^{logO/H) 



(All) 



where we define the gradient = A radial gradient in therefore naturally results in a non-zero metallicity 
gradient. It is instructive to examine the resulting evolution of metallicity gradients with time. Noting that the gas 
fraction was found to be approximately constant with radius in the lensed galaxies f iJ2.5l) . we make the simplifying 

assumption that ^^^fe^ = 0. In this case the rate of change of a metallicity gradient is given by 



dAz 
dt 



|A,A(logO/i/) 



= Afl 

= Ar 



jji^ogO/H) 



dA 



R 



dtdf',' ° ' ' dt df, 

1 (1 — a;) exp (x) — 1 



-{logo I H) 
cIAr 



(A12) 



1 -|- a; — exp (a;) 



ln(10)Mgas [exp(a;)-l]2 



dt ln(10)(l + /^)[exp(x) 



1 



A.l. Validity of the Model 

We now discuss important assumptions of the chemical evolution model, their validity, and quantify the implications 
for interpreting data using this model. There are several assumptions with which we are concerned. 

Perhaps the most important assumption is that total gas mass is constant fEquation lA2[) . Given that gas accretion 
and outflows are undoubtedly present, we expect the gas mass to fluctuate with time. However, we also expect the gas 
reservoir to be in equilibrium between accretion, star formation rate, and outflows. Cosmological simulations which 
incorporate outflows indeed show that the observed mass-metallic ity relation and its dependence on star formation 
rate can indeed be explained in terms of such an equilibrium (e.g. IDave et al.ll20lil) . In these simulations, departure 
from this equilibrium gives rise to scatter about the observed m ass-metallicity-SFR re lation, which is measured to be 
only 0.05 dex in metallicity among local star-forming galaxies ()Mannucci et al.ll2010D . Furthermore, if the gas mass 
remains constant, then the expected growth in stellar mass (see ^4.21) results in gas fractions /gas — 0.03 — 0.15 at 
z = for the lensed galaxies. This is in good agreement with direct measurements of gas fractions in local galaxies 
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with appropriate stellar masses logAf* ^ 10.7 (e.g. lYoung fc Scovillill991[ ). We therefore expect the assumption that 
Mg = to be approximately valid, resulting in an uncertainty of only ~ 0.05 dex in 12+log(0/H). 

Another critical assumption is that gas infall and outflow rates {fi, fo, fen) are a constant fraction of the SFR. This 
is clearly not physically true; these rates should therefore be interpreted as a time-averaged value over the formation 
history of the galaxy rather than instantaneous rates. Observationally it is well known that nearly all galaxies observed 
at high redshift drive outflows of gas, while outflows are much less common in their descendants at z ~ 0. Therefore fo 
must decrease with time, while /, and/or fem must decrease accordingly to maintain equilibrium. Therefore we expect 
that model-inferred fi, fo, and fen correspond to time- averaged values and are generally larger than the instantaneous 
rates. 

We have assumed in the model that accreted gas is divided into two components with metallicity Z — and Z — Zg. 
This is entirely for ease of interpretation in the context of this paper. In reality accreted gas will most likely have 
an average metallicity within the range < Z < Zg-, the values of fi and fem can be adjusted to reach the desired 
value. However for our analysis it is more convenient to think of the accreted gas as a pristine cosmological component 
and an enriched component, as formulated. If the enriched component arises from metal transport within a galaxy, 
then the assumption Z — Zg is equivalent to the condition that metallicity is a continuous function of position. This 
is necessarily true. The assumption that inf ailing gas has Z = 0, however, is not true. The inf ailing gas will have 
a metallicity Zi, introducing an additional term +ZiMi to the right side of Equation IA4I Functionally this can be 
though of as a modification to the yield, such that the effective yield is 

yeff = y + fiZi. (A13) 

The IGM oxygen abundance is measured to be [0/H] = —2.8 at z = 2.5 fSimc oe et al.l [200^ corresponding to 
Z, = 0.0015 Z© = 3 X 10-5. For yo = 0.0087, the effect of IGM metallicity is negligible (< 10% of yo) for /; < 30. For 
an infall rate fi = 7 corresponding to the highest acceptable value for the lensed galaxies (see Figure [TT|) , the effective 
yield differs from j/o by only 2.5%. For the redshifts and metallicities of interest in this paper, the assumption that 
accreted gas has Z = therefore introduces a negligible systematic error of < 2.5% or < 0.01 dex in metallicity. 

In summary, the assumptions used in our chemical evolution model are in reasonable agreement with the best 
current data and understanding of galaxy formation and of the IGM. The systematic uncertainties introduced by our 
assumptions are expected to be negligible (e.g. < 0.05 dex in metallicity). We caution that the inflow and outflow 
rates must be treated as time-averaged rather than instantaneous values. Values of fo and fi inferred from the model 
are expected to overestimate the current instantaneous rates. 
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